On the universality of the Discrete Nonlinear Schrodinger Equation 
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We address the universal applicability of the discrete nonlinear Schrodinger equation. By employ- 
ing an original but general top-down/bottom-up procedure based on symmetry analysis to the case 
of optical lattices, we derive the most widely applicable and the simplest possible model, revealing 
that the discrete nonlinear Schrodinger equation is "universally" fit to describe light propagation 
even in discrete tensorial nonlinear systems and in the presence of nonparaxial and vectorial effects. 
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I. INTRODUCTION 



The formulation of problems in Physics results essen- 
tially in theoretical models able to describe the prop- 
erties of a class of phenomena. Typically this consists 
of deriving a set of integro-differential equations, stem- 
ming from fundamental laws and shaving off all the as- 
sumptions which do not affect the predicted observables. 
The latter approach, known as Ockham's razor (OR) and 
summarized by the Latin sentence " entia non sunt multi- 
plicanda praeter necessitatem" (i.e. "entities should not 
be multiplicated beyond necessity"), is commonly em- 
ployed to discriminate between equally explanatory the- 
ories or competing systems of hypotheses. Furthermore 
the major goal of theories in physics should be the deriva- 
tion of "universal" models, based on the conjecture that 
various phenomena are ruled by common characteristic 
principles, thereby allowing one to apply the same de- 
scription to events taking place in different areas/fields. 
The question arising naturally is then: how can we be 
sure that a model really conforms to the Ockham's ra- 
zor or, equivalently, how can we guarantee that it en- 
tails both universal applicability and importance (see e.g. 
[H, 0,1 HQ)? "Universality" should be intended here as 
the model counterpart of "universality" in statistical me- 
chanics, the latter based on the concept that properties 
exist for and apply to a large class of systems indepen- 
dently of their dynamical details (see e.g. @, Chapter 21] 

and @, a a m m m m m m ) . 

Following the pioneering numerical experiments of Fermi, 
Pasta and Ulam on anharmonic lattices [HI [13], the 
study of energy dynamics in nonlinear discrete sys- 
tems has attracted growing attention in recent years. 
Thereby two major theories of the last half-century, 
namely the theory of integrable systems and the theory 
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of chaos, have been specialized to account for discrete- 
ness pil Ojl, |H [2l[. In this framework the noninte- 
grable discrete nonlinear Schrodinger equation (DNLS) 
has received significant interest [H [H, and references 
therein]. The DNLS is frequently encountered in several 
areas of science, including biology 24] , solid state physics 
[25l ] . Bose- Einstein condensates [26] and nonlinear op- 
tics [27| , the latter being particularly relevant in view of 
the mature techn olog y available for optical waveguides. 
[H, [H, Oil [m|, [sHHf. Despite its relevance the "univer- 
sality" of this model has never been addressed: in non- 
linear optics for instance the DNLS is derived assuming a 
paraxial regime, scalar Kerr nonlincarities and first order 
perturbations (see e.g. [33. C hapter 13], [3(| Chapter 11], 
[3lT pages 269-290], and [27|,[33|]). Such hypotheses left 
a number of open issues concerning large nonlinear re- 
sponses, large index perturbations (e.g., in deep gratings 
and in photonic crystals vectorial and nonparax- 

ial effects. The last, a subj ect of considerable interest 
in continuous media [H, [U [H, HtJ E(| , is of particular 
relevance in nonlinear optics owing to the recent progress 
in nanotechnology and the realizability of optical struc- 
tures with subwavelength features. Before proceeding, 
however, we want to state clearly an operative definition 
of universality for physical models. By "universality" 
we mean the widest applicability of a model in describing 
the largest family of phenomena and the various aspects 
of the dynamics they encompass. We aim hereby at veri- 
fying the universality of the DNLS equation in nonlinear 
optics, adopting an approach of general applicability. We 
employ the language of symmetries and Lie transforma- 
tion groups. They provide a powerful approach to the 
analysis , as witnessed by the large number of examples 
ranging from pure mathematics to chemistry, biology, op- 
tics, thermodynamics, solid state physics, quantum me- 
chanics and robotics |M El El El Ell E(| . In partic- 
ular symmetry methods and Lie transformation groups 
— originally developed by Sophus Lie in the 19th century 
EtJ — provide a systematic avenue to study nonlinear dif- 
ferential equations [4l|,|43|. Our analysis develops along a 
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two-fold top- down/bottom-up approach: in the first stage 
we perform a symmetry reduction of the problem stem- 
ming from fundamental laws and in the most widely ap- 
plicable case, i. e. full vectorial Maxwell's equations in 
the nonparaxial regime and in the presence of a tenso- 
rial nonlinearity. In the second phase we proceed via 
Lie symmetrie to solve a classification problem and to 
assemble the simplest equation sharing the properties of 
the model derived in the previous step. In comparison 
with existing literature [23, HH, [33| we remark that 
our approach deals with the large class of cases including 
nonparaxial vectorial propagation as well as a tensorial 
nonlinear response. This two-fold procedure yields the 
most widely applicable model with the simplest (i.e. the 
most elementary) structure possible: the latter is ensured 
by the classification (bottom-up analysis) and the former 
by the reduction step (top-down analysis). Furthermore 
it satisfies two fundamental requirements of predictive 
sciences: Ockham's razor and the principle of noncon- 
tradiction, the latter stated mathematically. We reveal 
that the DNLS is far more important than known to date 
because it is the simplest "universal" model of discrete 
nonlinear energy propagation, even in the presence of 
nonparaxial and vectorial effects, e.g. in optical lattices. 



with ([2]), yield Maxwell's equations: 

vxh =! d ' 
vxe =^ b < 



V-b 
Vd 



0, 
0. 



(4) 

(5) 

(G) 
(7) 



with d = do + p . In a one-dimensional lattice [27j the 
dielectric tensor is periodically modulated (with period 
A) along one direction, say y, and can be expanded as: 



e 



V = Y2^(x,y - nA) = V(r) + Ae(r) (8) 



with e (r) defining the canonical structure and Ae(r) an 
arbitrary linear perturbation. We then Fourier transform 
the potential: 



0(r;t)' 
a(r;i) 



duj 



$(r;w) 
A(r;u;) 



(9) 



and define the total Lagrangian (in the frequency do- 
main) La'- 



II. SYMMETRY REDUCTION ANALYSIS 

We begin by modeling the propagation of electromag- 
netic waves in a nonlinear medium using the action inte- 
gral T = J Cdrdt and the Lagrangian: 



£ = (d 



))e - h ■ b 



(1) 



with fields (e, d , h, b) defined by potentials (4>, a): 



-V0- 



d 



V x a = (2) 



the nonlinear polarization p NL describing a generic non- 
istantaneous tensorial Kerr perturbation of components 

pf L = X% ® eJ ® e?C ® e *> U.M G [1,2,3] = 
[x, y, z] (assuming Einstein summation over repeated in- 
dices) with x ( r ]t)ijkl the nonlinear susceptibility and 
<g) a convolution operator such that / <8> a <S> b ® c = 
US #id&dk/(r; *-6)a(r, &)&(r, &)c(r, 6) 

p8j . The Euler-Lagrange (EL) equations are derived as 
conservation laws through Noether's theorem [4j| by ex- 
ploiting the variational symmetries generated by the ba- 
sis {vi = <9/<90, v 2 = J2j d / da j} (J G [x,y,z]) [see e.g. 
El . Chapter 4, page 276]. We obtain: 



dC 

duL 



0, 



{a e [1,4]) (3) 



with the four-element potential u 
d/dj and u J a 



djU a . 



b,a x ,a y ,a z ), dj = 
The EL equations ([S]), together 



£ + £, 



C = D E* HB* 

£= (Dq + P NL /2)E* — HB 



(10) 



where C and C are the Lagrangian fields modeling the lin- 
ear canonical structure and the whole nonlinear medium, 
respectively; the former contains: 



D = e E, 

E = -V<l - iloA, 

and the latter: 

D = VE, 

E = -V3> - iujA, 

P NL = ^E 



B = 
B = V x A 



B = VH, 
B = V x A, 



(11) 



(12) 



with Setj = e Xini(r;w)^*£j(% + l/2)/2 [48]. From the 
group of rotations v = Ajd/dA j - A^djdA 1 + <l<9/<9$ - 
Qd/dQ + c.c, through Noether's symmetries, we obtain 
the conservation law: 

V(E* x H + E x H*) + tuj(e ( Ae + P NL )E* = (13) 

which can be regarded as a generalization of the Lorentz 
reciprocity theorem. If we apply the same procedure to 

the canonical structure only (i.e., for Ae = P NL = 0), we 
obtain: 



V(E* xH' + E'xH*) = 



(14) 



where (E, H) and (E', H') are two sets of solutions. The 
conservation laws (fT3")) - (fl"4")) support the expansion of the 
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fields (E, H) into an orthonormal eigenbasis provided by 
the canonical structure, resulting in a series of first-order 
differential equations completely equivalent to Maxwell's 
equations. In particular Eq. (jT5J) yields the orthogonal- 
ity between eigenvectors while Eq. (p~3|) allows us to cal- 
culate the expansion coefficients via orthogonality. The 
canonical eigensolutions are of the form: 



exp(i/3 v z) (15) 



"E„~ 








E vz 




. H "_ 


"( 




+ z 


Hvz 


) 



with subscript t denoting the transverse component, [3 V 
being the propagation eigenvalue, the spectrum of which 
encompasses a discrete real set (bound states or guided 
modes), a continuous real set (unbound states or ra- 
diation modes) and a continuous purely imaginary set 
(evanescent states) [50]. The orthogonality between the 
eigenstates is expressed by the density conservation of 
(fT4|) along z: 



dxdy{~E* vt x 



1/3, 



(16) 



Through (jTSJ) , the vector field (E,H) can be expanded 
by employing the transverse portion of the canonical 
eigenvectors, the longitudinal components following from 
Maxwell's equations. To perform this expansion, we be- 
gin by writing the transverse portions of the electromag- 
netic fields as: 



E t 
H, 



dvC v (z) 



Ei/t 



(17) 



with envelopes C v varying in the direction of propagation 
z. We then write the Euler-Lagrange equation of motions 
corresponding to the Lagrangian coordinate A z : 

V t xH( = iluco [(4 + Se zz )E z + Se zt • E t ] = (18) 

being 5e z t a vector of components 8e zt — [8e zx , Se zy ] and 
having expanded l/(e' + 5e zz ) = (l-8e zz /e')je' ' +0{8e 2 zz ) 
to first order in 8e zz (since Se zz <C e' ,36]). We find 
the expression of E z by direcly solving (fTgj) , keeping into 
account the transverse modal expansion (|17[) : 



dvC u {z) 



E vz i 



1 



K zz 
e' 



x e 



(19) 



being E vz — Vt x H . vt . We then write the z— component 
ofEq. ©: 



V( x E ( = -iujfiH z , 



(20) 



and solve for H z , keeping into account Eqs. (jTTJ), thus 
obtaining: 



Hz 



dvC v (z)H vz 



(21) 



being H uz = » Vt ^, / f ■ Summing up, our expansion in 



transverse modes is: 

E t * 
H, 



dvC v (z) 



E z 
E,, 



dvC v (z) 



E vz 



^(1 



E yi 



(22) 



(23) 



The integrals (fTT)) - (f2"3")) should be calculated over the 
whole spectrum of (3; however, since the power P carried 
by each eigenstate is P = T vv j\ we are interested 
here in the evolution of the bound states (i.e. modes 
with P G 3? and well confined within the canonical struc- 
ture). The latter evolution is obtained as follows. We 
first integrate Eq. (fT3|) across the plane (x,y): 



dxdy 



8 



(E* x H + E x H*) 



dxdy 



zw(e Ae + P 1NL )E 



0. (24) 



We then substitute Eqs. (|2"2"] ) -([2"5 )) into f2"i]). The contri- 
bution of the first integral of (|2"4"|) , keeping into account 
the normalization condition (1161). is: 



dxdy(E* tn x H tm + E tn x H* m ) 



7-Cm + ?>(ftm ~ Pn) 

az 



= -8~z C - ^ 



where we employed the indices n and m to indicate the 
canonical eigenvectors of the system. The substitution of 
Eqs. (j2"2"|) - ([2"3")) into the second integral of (j2"4"|) , after some 
straighfroward algebra, yields to the DNLS equation: 

d 

i-^^n + V'ra+l + lpn-1 + \lpn\ V>n = 0, (26) 

where we introduced the dimcnsionless £ = zK n _ n+ i and: 

K £ 



2K, 



n,n+l 



K 



n,n+l 



(27) 



with: 

K n ,n+1 



uje 



n, Z J - J n+l, Z 



dxdy 
x I AeEji^En+i^ 

= 5e E„ + zAe 

8e zx Ein,x + 8e zv E n ,y — £Se zz E n , z /e' 



(28) 
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where <jey = <J l o 3 (b l j + 1/2), [cr x ,a y ,a z ] 



[E n>x , E n>y ,eE n>z / e'] and n + h denotes the bound state 
of the canonical structure centered in y = hA. It is worth 
stressing that the reduction procedure based upon sym- 
metries of Maxwell's electrodynamics is exact; the only 
assumption used in deriving Eq. (|2"6"j) . in spite of the OR, 
is that guided modes are well confined. The latter hypoth- 
esis, stating that (3 n ^ (3 m Vm, n JEfy , allows us to neglect 
mismatched terms given by bound/unbound/evanescent 
states proportional to exp[z(/3„ — f3 m )z\ and guided modes 
with n + h > n + 2. The latter turns possible due to 
the exponential dependence of Kmn on the the distance 
between sites n and m [see e.g. Hi], Chapt. 13, page 522, 
Eq. 13.8-9], yielding an appreciable contribution only for 
nearest- neighbor channels n and nil. 
This analysis demonstrates that the discrete nonlinear 
Schrodinger equation models the dynamics in the gen- 
eral case of discrete nonparaxial vector propagation in 
lattices with a nonlinear tensorial response. 



III. CLASSIFICATION ANALYSIS 

The properties of the DNLS are established by its 
structure. In particular 
i) The two invariants of motion: 

n 

E = J2 + C^n+l + ^n| 4 (29) 

n 

where W physically corresponds to the power and H 
to the Hamiltonian. These quantities, generated by the 
symmetries {vi = J2 n iip n d/dip n + c.c.,v 2 = d/dz} of 
the Lagrangian C = 7^ hM*-^r + c.c.) — H , charac- 

as well as chaos 



terize nonlinear waves [33J as well as chaos [52 
ii) The existence of linear plane-wave solutions ip n — 
exp(if3z — iqn) with a (real) periodic dispersion relation 
/3 = 2 cos q [3l| • 

To obtain the simplest differential equation with the sym- 
metries expressed by {vi,V2} and compatible with ii) 
firstly we have to solve a classification problem, which 
yields all the equations with given symmetries. To this 
extent we generalize the theory developed in [Bli ] for 
difference schemes. Specifically, we start by consider- 
ing a general system of differential-difference equations 
Q = [Q^Q^Q 3 ] = involving nearest-neighbor (dis- 
crete) interactions: 



'Q 1 ' 




-dip- 
dri 




/ 2 (?7,x,gL 






<M, [ip-, 


<M) 


Q 2 




d± 
dr} 

d+ 








<t>+], [i>-, 


0-]) 


Q 3 






f 3 (v,x,d- 






<t>+], bP-, 


0-]) 



(30) 

with unknown functions /' (i € [1, 3]) depending up 
on the discrete variables [ip(x n , rf), <fi(x n , rj)] — [ip,<fi] = 
[i>,ip*] and their neighbors [i/>(x n ±i,rf),<j>(X n ±i,ri)] = 



[4 ! ±,4>±] = [V'i j V'dh] j defined in discrete [x, d+,c?_] = 
[x„, x„+i — x n , x n — x„_i] and continuous r\ = i£ spaces 
of independent variables. The first two of (|30|) describe 
the evolution of ip and <f> and the third defines the lat- 
tice. The introduction of the complex space r/ guarantees 
that tpandcp be real and the classification problem to be 
solved in the real domain 5i. To perform the classifi- 
cation, we need to find the structure of /' supporting 
the symmetries generated by the basis {vi = d/dr], v 2 = 
ipd/dip — cj)d/d4>} and satisfying a real dispersion relation. 
Each of the symmetries possessed by ([3TJ1) should satisfy 
the invariance criterion: 



=n=0, 



? J =o 



i 6 [1,2] 
J 6 [1,2,3] 



(31) 



with pr 1 the prolongation operator (of order 1) defined 
as: 



pr v = v + C(x+,y+) 
d 



d 
dx + 
d 



for a generic vector field: 

v = C— 

dx 

with: 

7 = [tp,4>}, 

7± = [i>±,<f>± 



+ ((x-,y-) 
d 



d 
dx- 



drj 



d 

'fry* 



+ cr l (x+,y+) 
(32) 

(33) 



d( dq 

drj drj dx 

drj 



dr d^j 
drj drj 



(34) 



The substitution of (f3"0"|) and (|3"2"|) into (|3 1 [) results in a 
system of equations which can be solved on the charac- 
teristics to yield the functional form of /' (i £ [1, 3]): 



-dip- 




dv 




d± 




dr] 




U+_ 





ipf x {x,d-,ili<l>, ip(f>±,H>±, ^f) 
-cj)f 2 (x,d-,ip(t), ip4>±,(jnp±, ■ 
d- f 3 (x,d_, tp(j), ip(f>± , <jyip± , ■ 



i± ij)± ■ 



(35) 



The simplest differential-difference equation encom- 
passed by ([33)) is obtained for f 3 — 1 (i.e., for d + = 
d- = d) and for / , f 2 depending on just one nonlinear 
product ip<p and one couple of linear terms ty±/ip, 4>±/'4'] 
, the latter giving a real dispersion relation as required 
by ii): 



dip 
drj 



4>+ +ip- + 4> 2 4> 



i>+ + 4>- + (p 2 ip 



drj 



(36) 



The factor 1 jd guarantees the existence of the continuous 
limit of (|3"6"|) for d — > 0, [ip(r]), <f>(r))] = d[u(x,r)),v(x,r))] 
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and x = nd. However, in our dimensionless example we 
can always set d = 1. Finally, going back to the original 
variables [£, tp, ip*], we obtain the DNLS equation ([26]) . 
The solution of the classification problem therefore 
demonstrates that the discrete nonlinear Schrodinger 
equation has the simplest structure compatible with its 
properties. 

In summary our top-down/bottom- up analysis shows 
that: 

I) The DNLS is "universal" because it is able to model 
discrete nonparaxial vector propagation in lattices with 
a tensorial nonlincarity. This stems from fundamental 
principles and exploits the conservation laws arising 
from Noether's symmetries [Eqs. (fT3| -([T4 |) ]. 

II) The DNLS is "simple" because it contains the 
smallest number of terms that provide a functional 
form compatible with its properties. This stems from 
solving the classification problem [Eqs. [35] by means of 
a Lie-symmetry analysis. 

We can therefore state that the DNLS equation is a 
general model, the importance of which goes far beyond 
the limits within which it is usually derived, being 
the simplest model able to describe nonparaxial vector 
energy propagation in a discrete medium with a tensorial 
nonlinear response. It is worth underlining that this 
result was obtained by simply exploiting symmetries, 
i.e. a general approach that can be extended to all the 
physical contexts, where yet model "universality" has 
not been assessed. To provide an example we consider 
the case of two (or more) dimensional optical lattices 
(see e.g. \Kk. \Eb\. \Edj). Our theory, although applied in 



the present context to the case of one dimensional optical 
lattices, can be straightforwardly generalized to two (or 
even more) dimensional systems, allowing to investigate 
the universality of higher dimensional physical models. 



IV. CONCLUSIONS 

In conclusion we undertook a general procedure to es- 
tablish both the significance and the "universality" of 
equations describing physical events, attributing such 
features to those which model the most general dynam- 
ics in the most elementary form possible. We demon- 
strated our approach in the case of optical lattices, re- 
vealing that the discrete nonlinear Schrodinger equation, 
previously considered in the specific context of scalar, 
paraxial propagation in weakly modulated media, is far 
more "universal" and important as it describes discrete 
nonparaxial vector energy propagation in materials with 
a tensorial Kerr response. The general character of the 
outlined procedure as well as the contemporary interest 
of discrete nonlinear systems are expected to contribute 
towards a better comprehension of model validity and 
applicability in physics. 
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